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Abstract 

Electron lyatems which have low dimensional properties have been constructed by squeez- 
ing the motion in zero, one or two-direction. An isolated quantum dot is modelled by a 
potential box with delta-profiled, penetrable potential walls embedded in a large outer box 
with infinitely high potential walls which represent the work function with respect to vac- 
uum. We show the smooth crossover of the density of states from the three-dimension to 
quasi-zero-dimensional electron gas. 


1 Introduction 

Quantum wires and quantum boxes with three-dimensionally confined electrons constitute a con- 
siderable part of recent semiconductor research ( 1 , 2]. To study the optical properties of these 
systems, one should investigate the density of states (DOS) carefully, because the change in the 
density of states affects directly the optical properties of these structures as a result of reduced 
dimensionality. 

The DOS of a low dimensional electron gas(LDEG) in the presence of magnetic field has 
been discussed in many literatures measuring the magnetocapacitance [3, 4]. Furthermore an 
electrical confinement which is usually controlled by (alternate) gate voltage [5] and, so called 
the illumination method [6] are used to get a LDEG. The etched silicon filaments also discussed 
recently as quantum wires or quantum dots [7]. But the DOS of a LDEG of confined electrons 
in small space which is constructed by reducing the size of the confinement is not discussed 
frequently, see ref. (1). A typical example of an ideal system having QOD charactor is that of 
electron confined in a quantum box with impenetrable potential barriers. Despite of the large 
number of studies on quantum wire and quantum box structures up to date, we have not found 
research on the crossover of the DOS from a three dimensional DOS to a quasi- zero dimensional 
DOS. In Section 2, to illustrate the formation of a quasi-one-dimensional electron gas(QlDEG) 
using the classical electrostatic method, a simple metal-insulator-semiconductor(MIS) structure 
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with very many parallel gate electrodes has been treated by making use of the conformal mapping 
method In section 3, we consider a rather artificial quantum box structure, so called the three 
directional double-barrier resonant-tunneling structures(DBRTS) to study quantum mechanically, 
and have calculated the local density of states and the global density of states. In section 4, the 
crossovers of the DOS is calculated. Especially, we reveal the crossovers of the DOS from 3D to 
QOD. 

2 Construction of the very many parallel quantum wires 

To study electrical properties of a quantum wire, we first start with a quasi-two-dimensional 
electron gas(Q2DEG) at simple metal-insulator- semiconductor(MIS) structure. A Q2DEG system 
with many parallel gate electrodes is shown in Fig-1 in which an electron gas is confined to the 
x-z plane. We actually try to confine the electrons in the z-direction as well to form a QlDEG 
system. 
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Fig.l A Structure of symetric gate arrays. Fig.2 The boundary condition in w- 
plane. 

To calculate the charge(density) distribution at the MIS interface (y=0 plane) to see the 
formation of a QlDEG, we will use the conformal mapping method which is useful especially ior 
two-dimensional problems and we assume that significant changes in the electrode potential (and 
thus is density in the channel) cause only a slight change in the near junction band bending. 
This type of approximation has been used by Shik [8] to calculate various properties of the MIS 

structure. 

The problem is solving the Laplace equation in insulator region. 


d 2 V 
dx 2 + dy 2 


( 1 ) 


with the boundary condition; 
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At y « d, * |y _4 it alternate gate potential Vi and V 7 and at y = 0, * |y. 0 « constant, i.e. 
equipotential surface. The next step is getting the distribution of the surface carrier density. 

B#( * } * Axtdy lr "° ’ (2) 

We solve the problem by taking the following conformal transformation (9j. 

W * U + «'V = e wtfi ,z = x + iy 

where d is the thickness of the insulator, Now the insulator region is mapped on the upper half 
plane and the boundary condition is given as in Fig.2. 

After getting the potential which satisfy the boundary condition, we now get the electron 
density distribution n(*) analytically from Eq.(2): 

n(x) - ^[1 + M/V,)- 1 B-ir 1 < sinh((2n - l)xa/2d) 

* We nml 

[cosh (tx/d) + cosh((2n - l)/2</)] >] (3) 

where a is the gate interval and nt is the number of gates. 



Fig.3 The surface charge density vs. position. 

Fig.4 Crossover of the global DOS from 3D to 2D in the range U\ — 0 to U 1 — 20, 
as a function of E/Eo. Here Ui takes the values 0, 2, 12, 16, 20. 

A typical density distribution is illustrated in Fig.3 where we can see immediately the many 
parallel QlDEG (eventhough we show here only two wires). The one dimensional electron density 
of the order of 10 6 /cm is obtained for the typical operating gate voltage when a is 1000 A. We 
also investigated the case of anti-symmetric gate voltage. Similar results have been obtained but 
in symmetric case it is easy to construct one dimensional electron channels especially for smaller 
number of gates. In our calculation we took a/d * 1, which satisfy the first approximation, y = 0 
plane is equipotential. 
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3 The quantum box model and the DOS 

Now we come back to quantum system with a delta-profiled quantum box. Usually a quantum dot 
is an dement of the array of quantum dots. But the interaction among quantum dots decreases 
rapidly with increasing dot separation [10] and is unimpotant for the usual experimental situation 
[11]. 

Therefore a separated single quantum box is taken for our study. For a rectangular wire(Z, v(r) < 
Ac), Arora and others[12] used the impenetrable potential walls, but in this paper we consider 
three sets of penetrable barriers. We start our calculations with the model, i.e., the typical three 
directional DBRTS, which consists of two thin(~ 50A)Al,Gai_,Aj layers, separated by a thin 
GaAs layer along all three directions. The potential is expressed by 

V(x,y, z) ■ V x $(x + a) + 6{x - a) + V>6(y + 5) + 6(y -b) + V£d(z + c) + 6(z - c). (4) 

In this potential, the six Al t Ga\- f Aa potential barriers have been replaced by ^-functions with 
strengthes Vi, Vj and Vs in the x,y, and x direction, respectively. The parameters V(t = 1,2,3) 
are given by 

Vi « *AV a (5) 

where di are the barrier widths and A Vi, are the conduction-band discontinuities. In order to deal 
with finite density of states,[13], we must place our structure in a large impenetrable rigid box 
extending from —L/2 to L/ 2. With proper boundary conditions[14], the Schrodinger equation is 
separable, we can write the wave function in the product form 

= = ( 6 ) 

The separated wave functions, satisfy the reduced equations. 

+ pme/ft’p, - Vi]#. = 0 (7) 

with 

E = t,E, ( 8 ) 

isl 

Here E is the total energy corresponding to the Hamiltonian H and Ei(i = x, y, z) is the energy 
eigenvalue of 

The local density of states in the DBRTS has been obtained in various cases [15]. It is defined 
as a function of r * (x, y, x) and E by 

N(x,y,x’,E) - -(2/»)/mG(r,r' ; E) 

- 2££EE I ♦.*.(*) 1*1 *Wv) l a l **.(*) I J 6(E - £*), (9) 

ajh ** *» 

where the factor of 2 implies spin degeneracy, G{t^\E) is the single particle Green's function, 
and a,0, and 7 (= e or 0) label state parity. Next we consider the global DOS N(E). It can be 
calculated by taking the integration over the box volume, 

N(E) « 8 /* dx jf* dy J‘ dzN(x, y, z; E). ( 10 1 
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The Amplitude of the wave function inside the well for both even and odd parities of x, y, and z 
componenett are given elsewhere[14]. N(E) can be rewritten as follows: 

N{E) * (2/s” 3 ) jT°dpi[G,(pi) + G*(pi) + (G.Oh) - <7,(p 1 ))sin(2pi)/2pi] 

* dPilG»(p>) + G»{Pi) + - G*(l>a))«i“(2p*)/2pj] 

* dp i [G,(pi) + G 0 {pa) + {G,ipi)-G,(p i ))mii( 2 p a )/ 2 p 3 ]S(E - E k ). (11) 

The properties of functions G,(p,) and G 0 (j>i) are already revealed in Ref.14. 

When we take appropriate limiting cases, the Eq.(ll) recovers all the well-known expressions 
of the DOS of 3D, 2D, ID, and OD. since the calculations are straightforward, we haven’t repeated 
here. 

4 Crossovers of the density of states 

Now we consider crossovers of the global DOS from a high dimension to a low dimension. 

4.1 From 3D to 2D 

This case may happen when two of three potentials U{»( see the reference 15) approach zero, while 
the remainder varies from zero, i.e., 3D case, to infinity, i.e., 2D case. The Eq.(l 1) can be modified 

as 

MAM + <?.M + (G.M - <J.(p,)).in(2p 1 )/2p l ] (12) 

The result of the numerical behavior of Eq.(12) is shown in Fig.4 and indicates the transition of 
the DOS from 3D to 2D. In this case we take U t * * 0, Uj changes from 0 to 20, and E/E 0 

varies from 0 to 8. Higher values of U\ correspond to a staircase-like 2D behavior which shows 
steps at E/E a * P with / * 1, 2, 3, • • •. 

4.2 From 2D to ID 

This corresponds to the case of f/j going to zero, U\ to infinity, and U% varying from zero, i.e., 2D 
case, to infinity, ID case So G(ps) » 1 ,<?.( r ) « »£*(pi -(/ + l/2)w),G.(pi) = »£*(pi - (/ + 
1)t). Then we can get the modified equation of N(E) as follows ; 

» J /l l ~ f/HB-B,) 

j—N(E) - £/ 4k[G.{?i) + GM + {GM-GM)*, in(2 ft )/2 ft ] 

SCm e mmO J0 

x [1/[(jt/2 a)’E/E. - ((m + 1/2 )x/a) J - (p a /5) , ) 1/1 
+ l/[(*/2a) 2 E/E 0 - (m + l)x/a J - (p,/*) 3 ] 1 ' 8 ]. (13) 

The Fig.5 shows the graphical result, that is, the crossover of the global DOS from 2D to ID. In 
this case, for the sake of convenience, we take t s (,[/, « 0 and U x — * oo, Uj takes the values 
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of 0, 2, 8, 16, 20, and E/E, is taken from 0.0 to 8.0. One can see that the steps at E/E, = 1, 
4, 9 ] •••, that is, the two dimensional band edges, are shifted to peaks at E/E, » 2, 5, 8, 10, 

• • •, i.e., the one dimensional band edges. As U 3 increases, the motion of confined electrons along 
the* y axis atarts to s hr ink and is guided only along the * axis. This weak additional confinement 
shifts the 2D band edges towards higher energies and finally the typical ID characteristics of the 
DOS comes to bevisualized. Higher values of t/j correspond to increased sharp peaks of the DOS 
of the 1-D quantum wire case, which are in good agreement with those of Arakawa and Sakaki 
[13] and of Tsang [16]. The values at E/E, * 5 and 10 are roughly twice those at E/ E, = 2 and 
8, respectively, which comes from the double degeneracy of the eigenstates. Similar discussions 
were treated by Berggrent and Newson [17] in the case of the 2D electrons in the presence of a 
magnetic field. 



Fig.5 Crossover of the global DOS from 2D to ID. Here we take U\ ooandVj — 
0, 2, 8, 16, 20. Higher values of U% correspond to a sawtoothlike ID behavior. 

Fig.6 Crossover of the global DOS from ID to 0D. Here we take U\ — Uj ~ 00 
and U 3 * 0, 2, 10, 20, 80. Higher values of U 3 correspond to a sharp line shape 0D 
behavior. 


4.3 From ID to OD 


In this case, we take both U t and U 3 to be infinity, and U 3 vary from zero, ID case, to infinity, 
0D case. Then Eq.(ll) becomes 


N(E) 


OO 

<»/*> E £ 

l.mml ,/0 


ME/S.) 1 '* 


dptlGM + G.(ps) + (G.(ps) - Go(pj)) Sin(2p3)/2p3] 


x 6[E - (fi 3 /2m e )(/T/2o) a + (mx/24) J + (ps/c) J 

A C.(0 4- g.(t) + (g.ft) - i 

(2 E.)N(E) = E [E/E,-P-Tn*Y n 


(14) 

(15) 
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where t = (*/2)(E/E 0 — P — 

Fig.6 shows the transition of the global DOS from ID to OD. For the sake of convenience, we 
put a m b * c. The Eq.(15) recovers the well-known DOS of a quantum dot [16], when we take 
U 3 to be infinity. Sawtooth type maximums at E/E 0 = 2, 5, 8, 10, • • •, are now moved to the 
positions at E/E 0 * 3, 6, 9, 11, • ■ as the strength of f/j increases. When the confining potential 
increase, both primary peaks and secondary ones appear, which reflect the coexistance of ID 
and OD behavior. The secondary peak with the lower energy is a reminescence of the ID DOS 
shifted towards higher energy due to the additional confinement, and the pri mar y peak (higher 
energetic peak) arises from quasi OD states. Bacause the differences between peaks are so high, we 
used different scales for the DOS axis ranging from 1.0 to over 2000. The DOS clearly shows the 
potential strength (U,) dependence of the spatial quantization through E 0 . This kind of secondary 
peaks are also shown in many experimental data of a transport measurement [18]. We know that 
the electron systems used in above experiments are in an intermediate state between ID and 0D, 
because the potential strengths are not infinitely high. 

We believe that this kind of DOS transition which shows intermediate states will also occur in 
real systems where, for example, the barriers have finite widths. For barriers with finite thickness, 
the effective mass of the electron changes in passing from the quantum-well region (GaAs) to the 
barrier regions (AlGaAs) of the structure. Bruno and Bahder [15] have considered this for the one 
directional DBRTS case and showed that the DOS at the low-energy subband edges is higher than 
the DOS at the same energies in the absences of barriers (for delta-profiled barriers). In our case, 
we can estimate that our result for the DOS will be increased a bit upward at the same energies 
because of the additive form of the potential which we have taken. 

In this paper, first we showed as an example the formation of many electron wires using the 
conformal mapping method. Next, considering a penetrable quantum box, with a volume of a 
axfocc, in a very large rigid box of volume L 3 , we calculated the general form of the local and 
global DOS. 

The merit of this model is as follows : 

1) the model is simple to handle and easy to calculate analytically, 

2) in this model, one can recover the results of all the limiting cases of the 3D, 2D, ID, and 
°D, 

3) starting from one equation we can discuss all three crossover cases. 
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